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Abstract 

This paper deals with a new class of random nights X d (i), t > 0, defined in the real space 
M. d ,d > 2, characterized by non-uniform probability distributions on the multidimensional 
sphere. These random motions differ from similar models appeared in literature which take 
directions according to the uniform law. The family of angular probability distributions 
introduced in this paper depends on a parameter v > which gives the level of drift of the 
motion. Furthermore, we assume that the number of changes of direction performed by the 
random flight is fixed. The time lengths between two consecutive changes of orientation 
have joint probability distribution given by a Dirichlet density function. 

The analysis of 'K d (t),t > 0, is not an easy task, because it involves the calculation of 
integrals which are not always solvable. Therefore, we analyze the random flight X^(i),i > 
0, obtained as projection onto the lower spaces E m , m < d, of the original random motion 
in R d . Then we get the probability distribution of X^(t),t > 0. 

Although, in its general framework, the analysis of X d (i),i > 0, is very complicated, for 
some values of u, we can provide some results on the process. Indeed, for v — 1, we obtain 
the characteristic function of the random flight moving in R d . Furthermore, by inverting 
the characteristic function, we are able to give the analytic form (up to some constants) of 
the probability distribution of X d (t),i > 0. 

Keywords: Bessel functions, Dirichlet distributions, hyperspherical coordinates, isotropic 
random motions, non-uniform distributions on the sphere. 

1 Introduction 

The random flights have been introduced for describing the real motions with finite speed. The 
original formulation of the random flight problem is due to Pearson, which considers a random 
walk with fixed and constant steps. Indeed, the Pearson's model deals with a random walker 
moving in the plane in straight lines with fixed length and turning through any angle whatever. 
The main object of interest is the position reached by the random walker after a fixed number 
of the steps. 

Over the years many researchers have proposed generalizations of the previous Pearson's walk 
randomizing the spatial displacements. In particular, the random flights have been analyzed in- 
dependently by several authors starting from the same two main assumptions. The first one 
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concerns the directions which are supposed uniformly distributed on the sphere. Furthermore, 
the length of the intervals between two consecutive changes of direction is an exponential ran- 
dom variable. Therefore, there exists an underlying homogeneous Poisson process governing the 
changes of direction. The reader can consult, for instance, Stadje (1987), Masoliver et al. (1993), 
Franceschetti (2007), Orsingher and De Gregorio (2007), Garcia-Pelayo (2008). A planar random 
flight with random time has also been studied in Beghin and Orsingher (2009). 

Exponential times are not the best choice for many important applications in physics, biol- 
ogy, and engineering, since they assign high probability mass to short intervals. For this reason, 
recently, the random flight problem has been tackled modifying the assumption on the exponen- 
tial intertimes. For example, Beghin and Orsingher (2010) introduced a random motion which 
changes direction at even- valued Poisson events. In other words this model assumes that the 
time between successive deviations is a Gamma random variable. It can also be interpreted 
as the motion of particles that can hazardously collide with obstacles of different size, some 
of which are capable of deviating the motion. Very recently, multidimensional random walks 
with Gamma intertimes have been also taken into account by Le Caer (2011) and Pogorui and 
Rodriguez-Dagnino (2011). Le Caer (2010) and De Gregorio and Orsingher (2011), considered 
the joint distribution of the time displacements as Dirichlet random variables with parameters 
depending on the space in which the walker performs its motion. By means of the Dirichlet law 
we are able to assign higher probability to time displacement with intermediate length. The 
Dirichlet density function in somehow generalizes the uniform law and permits to explicit for 
each space the exact distribution of the position reached by the motion at time t. 

In the cited papers, it is assumed that the directions are uniformly distributed on the sphere. 
This last assumption is not completely appropriate. Indeed, the real motions are persistent 
and tend to move along the same direction. Therefore, it seems not exactly realistic to have 
the same probability to move along each direction of the space when the motion chooses a new 
orientation. The aim of this paper is to analyze a spherically asymmetric random motion in 
R d , that is with non-uniformly distributed directions. Essentially, the literature is devoted to 
analyze random walks with uniform distribution, while on the random flights with non-uniform 
distributed angles few references exist: see for instance Grosjean (1953) and Barber (1993). For 
this reason we believe that this topic is interesting and merits further investigation. 

For the sake of completeness, we point out that other multidimensional random models with 
finite velocity, different from the random flights, have been proposed in literature: see, for 
example, Orsingher (2000), Samoilenko (2001), Di Crescenzo (2002), Lachal (2006) and Lachal 
et al. (2006). 

We describe the random flight model analyzed in this paper. Let us consider a particle which 
starts from the origin of R d , d > 2, and performs its motion with a constant velocity c > 0. 
Hereafter, we assume that in the time interval [0, f], a number n, with n > 1, of changes of 
direction of the motion is recorded. We suppose that the instants at which the random walker 
changes direction are < t\ < t2 < • • • < t n < t, with to = 0, t n +i = t, and denote the length 
of time separating these instants by tj = tj — tj-i, 1 < j < n + 1 where < Tj < t — X)fe=o T k, 
1 < j < n, and t„ +1 = t - Y!j=i T j- 

The directions of the particle are represented by the points on the surface of the ti-dimensional 
hypersphere with radius one. We assume that the random vector 6= (6\, ...,64-2, 4>), represent- 
ing the orientation of the particle, has density function given by 

<W,„O?i,02, ...,0 d -2,<P) = + m sin 2 '-+ rf - 2 6 1 sin 2 -+ d - 3 fl 2 ■ ■ ■ sin 2 ^ 1 9 d _ 2 sin 2 - (1.1) 

2n — T{v + 1/2) 

with 0j € [0, = 1, 2, ...,d— 2, <f> £ [0, 2tt] and v > 0. Hence, the particle chooses the direction 
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Figure 1: The picture on the left represents the behavior of density function g2,i{x) — - sin 2 x 
with x G [0, 2tt] . On the right it is depicted the density g$,i{x, y) = A sin 3 x sin 2 y with x G [0, 7r] 
and y € [0, 27r]. 



not uniformly on the sphere, but following an asymmetric density law. Furthermore, the particle 
chooses a new direction independently from the previous one. In Figure [T] we have depicted the 



behavior of (1.1) for v = 1 and d = 2,3. From the picture on the left we can observe that the 



function (|1.1[) is bimodal and assigns high probability mass near the points | and |7r. Similar 



considerations hold for the picture related to d = 3. Clearly when v increases the function g 
tends to concentrate around the points ^ and |7r. For v = 0we obtain the uniform distribution 
over the d-dimensional hypersphere, that is 



T(d/2) 
2ni 



sin^sin^ 3 • ■ • sin 0^-2- 



Another important assumption is the following: the vector r = (r l5 r 2 , t„) has distribution 
given by 

, , , r((n + l)(2^ + rf-l)) mt\r^ +d - 2 

Jd,v{Tl,...,T n ) (T(2 V + d- l)) n+1 i(n+l)(2»+d-X)-l ' 

which represents a rescaled Dirichlet random variable, with parameters (2v + d — 1, 1v + d — 
1), d > 2. 

Let us denote with X rf (t),£ > 0, the process representing the position reached, at time t, 
by the particle following the random rules above described. In the rest of the paper we are 
going to analyze the random flights defined by X d (t),t > 0, and their probabilistic character- 
istics represent the main object of our analysis. Recalling that the motion develops at con- 
stant velocity c, and exploiting the hyperspherical coordinates, the d-dimensional random flight 



Figure 2: A sample path of the three-dimensional random flight with its projection onto the 
plane. 



X rf (t) = (Xi(t), Xd(t)),t > 0, has components equal to 

n+1 

X d {t) = T fc sin6»i ife sin 6 2 ,k ■ ■ • sin0d-2,fc sin^) fe 
fc=i 

n+1 

-X"d-i (*) = c r fc sin 6»i ife sin Q 2 m • • • sin # d _ 2 ,fc cos </> fc 



k=l 

(1.3) 

n+1 

-^aCO = c^rfc sin (?i,fc cos g 2 ,fc 



fc=l 

n+1 



^l(0 = C^TfeCOsgi.fc. 



fc=l 



The probability distribution of random flight X rf (t),f > 0, is concentrated inside the d- 
dimcnsional hypersphere with radius ci, which we will indicate by H^ t — {x d : \\x d \\ < ct}, 



where x d = (xx,X2,—,%d) an d || • || represents the Euclidean norm. The density law (l.ll leads 
to a drift effect in the random motion. Indeed, since the directions are not chosen with the same 
probability the particle tends to move in some portions of the space R d with higher probability. 
Clearly the parameter v defines the degree of drift of the motion. Therefore, if the value of v 
increases, the particle will tend to move along the same part of the space. This behavior of the 
random motion seems to be consistent with the persistence of the real motions. 



The sample paths described by the moving particle having random position (1.3) appear as 
straight lines with sharp turns and each steps are randomly oriented and with random length 
(see Figure 2). 

The drawbacks emerging in the study of the random flights defined in this paper are discussed 
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in Section 2. Therefore instead of studying directly X d (i), t > 0, we analyze its projections onto 
the lower space R m , m < d. This approach leads to another random motions for which we obtain 
the probability distributions. 

For v = 1, we provide the characteristic function for X d (t),i > (Section 3). Furthermore, 
we infer a closed- form expression for the density function of the random flight with v = 1. For 
n = 1, 2, it is possible to explicit completely the above density function. 

The Appendix contains results on the integrals involving Bessel functions, which assume a 
crucial role in this paper. 



2 Motions related to the random flights 

In order to study the process X. d (t),t > 0, introduced in the previous Section, we take into 
account the same approach developed in Orsingher and De Gregorio (2007) and De Gregorio 
and Orsingher (2011). The first step consists in the calculation of the characteristic function of 
X d (i), t > 0, which plays a key role in our analysis. Let us indicate by < •, • > the scalar product, 
then we have that 



•t-n rt-ELi Tk 

dn I dr 2 --- d,T n f ( i,v(n,...,T n )Xn(Ti,...,T n ,a d ), (2.1) 



where a d = {a\, a d ) and 

^(t 17 ..., T„,a d ) (2.2) 

t>2n p2iT 



■>n+l 



/ dOis ■ I dQi, n+1 ■ I dQ d -2.\ ■ I d9 d - 2 ,n+i / #1 ■ • • / 
Jo Jo Jo Jo Jo Jo 

n+l ( 

Y\ l exp{icT k (a d sm9 ltk sm9 2 ,k ■ ■ • smd d ^ 2 ,k sm(f) k + ctd-i sin9 lyk s'm9 2 ,k ■ ■ • sin9 d ^ 2 ,k cos0 fc 



fe=i 



+ •• • + a 2 sin9 hk cos9 2tk + a x cos^)} ^}" + d ' 2 ' sin^' y + d - 2 ■ • • sin^+i . 

27r^r(i/ + 1/2) 

The integrals with respect to the n + l angles <p k are of the form 
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g iz(a cos <p+P sin <j>) 







which seems to have an explicit solution only for particular values of v (see Appendix) . Thus the 
choice of the non-uniform family of distributions (1.1 1 for the directions of the random motion, 
makes the analysis of the process, at least in its general setting, very complicated. Clearly for 
v = we get the uniform case studied in the paper mentioned in the Introduction. 

For this reason, instead of studying directly X<j(t), we deal with the random process X^f) = 
(Xi(t), X m (t)),t > 0, namely the projection of X rf (f) onto the lower space R m , with 1 < m < 
d, having components equal to 

n+l 

sin $i,fc sin $2,fe ' ' ' sin $m— cos 

fc=i 
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n+l 



X m -x(t) = c^T fc sin 9 ltk sin 6> 2 , k ■ • • sin0 m _ 2 ,fc cos0 m _i,fc 



fc=i 



(2.3) 



n+l 



^(*) = c Tfc sin g lt fc cos fl 2 ,, 
fc=i 

n+l 

= T k cos 6>i,fc. 



fc=l 



The vector (61,62, ...,6m) (with $d_i = cf>) appearing in (2.3), has distribution given by the 
marginal density of <?d,y (01, $2, @d-2, 4>), that is 



9d, v {Q-i 



P7T P 7T f2lT 

,6 2 ,--,6 m )= I d6 m+ f - I d6 d -2 I d(f>g 
Jo Jo Jo 
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d,v\Vl,V2, — ) Od-2, (P) 

m = d — 1 



r{v+d/2) sin9 2u+d-2 . .. S i n 02u+d- m -l < d _ 1 

7rTrr(i/+^=S) 1 m 



(2.4) 



Therefore, X^(i) can be interpreted as the shadow of ( |1.3[ ) in the lower spaces (see for instance 
Figure 2) . In other words, if we observe in R m the particle moving in M. d according to the random 



rules of X d (£), we perceive a random flight with components (2.3 1 having vector velocity 



/ c sin 0i sin 6*2 • ■ • sin6> m _i cos0 m \ 
c sin 61 sin 62 ■ ■ ■ sin0 m _ 2 cos0 m _i 



c sm t>i cos 6 2 
y c cos 0i 



/ 



with random intensity \v\. Then, in what follows we analyze the random motion X^(t),i > 0, 



function of (2.3) 



with orientations distributed according to (2.4 1. Our first result concerns the characteristic 



Theorem 1. The characteristic function of X m (t ) , t > 0, is equal to 



r, 



2 H±i (2l/+d _ 1) _i r(2 ±l (2 ^ + d _ + I) 

J n + l 



{ct\\a m \\)^^+ d -^- 1 2 

(-l) fc i x ^2k+n 



where J^(x) — J2k=o k\r(k+fi+i) \ 2 
Proof. The characteristic function of X^ (t ) , t > 0, becomes 



( 2l ,+d-l)-l( C *IKJI)> 

, n € K, is i/ie Bessel function. 



(2.5) 



^(a m ) = s{c*<a.^.W>} 



t /■*— Tl 

/ dr 2 ■ 

Jo 



dr n f d ,u{ri, ...,T n ) 
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with 



(n, ...,t„ 

, Q L7n 

) (2-6) 

/»7T f*7T f*7T 

ddi t i ■■■ d9i :n+1 ■ ■ ■ d6 mjl ■ ■ ■ I d9 m ^ n+ i 



n+l 



Y]_ < exp{icT k (a m sin0i,& sm6 2 ,k ■ ■ • sin# m _i )fc cos0 mi fc H \- a 2 sin9 lyk cos9 2 ,k + ui cosfli,*)} 



k=l 



T(v + d/2) 2u+d _ 2 ■ (fiv+d-m-1 



x — = ; sin ft, ■ ■ ■ smt t 



Since 



V2/ / O I2COS0 ■ 21/ D „ I ., I 



J » & = r(^ + i)r(i) Jo e " C ° S s ^ Re V + 2 ' > °' ! 2 7 ' 



we observe that, after the integrations with respect to 9 m _ k , k = 1, n + 1, (2.6 1 becomes 

/•7T /*7T /•TT 

• • • / d#l, n +l ' ■ ■ / d&m-l,l ' ' ' / ^#m-l,n+l 





n+l 



< cxp{icr fc (a m _i sm#^ fe cos9 2 ,k • • ■ sin# m _ 2 ,fc cos0 m _i,fc H h«2 sin#i jfe cos0 2 ,& + ai cosfli,*)} 

fc=i I 

sin ^- 2 --. sin e!Ji 



I— m 



2i/+d— m — 1 



TT f r(i/ + d/2) 

TT S ,„ =7 / exp{zcr fc (a m sin0i fe cos0 2 ,fc • • ■ sin0 TO _i fc cos0 ro , fe )}sm0, 

/•7T /"7T /*7T Z 1 ^ 

: / d9\s ' ' ' / d01,n+l ' ' ' / d9 m -l,l ■ ■ ■ d0 m _i.„+i 

Jo Jo Jo Jo 

n+l ( 

Y]_ < exp{icr fe (a m „i sin0i ife cos9 2 . k ■ ■ ■ sin0 m _ 2 ,fc cos 0^-1,* H h«2 sin0 1)fe cos0 2 , fe + a% cos0i ifc )} 



fe=i k 

r(i/ + d/2) . g+ri-2- "-y- 1 . v+d-m-i 

X 55=1 Sln ^l,fc •■• Sm ^m-l,fc 

TT 2 

/ 2 n ! 



x J (cr fc a m sin ^i.fc sin 6» 2 .fc • • • sin0 m _i ifc ) >. 

\cT k a rn J '+2 j 

We are able to perform all the (to — l)(n + 1) integrations with respect to the angles 9i. k , 1 < 
i < m—l,k = 1, n+ 1, by applying successively the formula (4.111 in the Appendix. The 
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integration with respect to 6 m —i t k, k = 1, n + 1 yields 

n+l 



Jfl TT „icT t a m _i sin6»i t—sin m _ 2 i cos# m _ik 

"fm-l,n+l XX e 

fc=l 

1 

d — m — 1 



-l.k 



I d-m-i (cr fc a m sin fli |fc sin g 2 ,fc • • ■ sin m -i,&) 

n+l 



n 

fe=i 



x — , d-m-i Jy+ (cTfeCtm singi.fc sin fl 2 fc ■ ■ • sin6» m _i, fc )d6l m -i,fc I 

n+l r „7r/2 

= II 1 2 / cos(cT fe a m _i sin 6»i. fe • • • sm9 m - 2 ,k cos6> m _i jfe ) 

sin6l m-i,fe t , , d-m-i ^ + d-m-i (cr fe Q; m singi fcSinfl? fc • ■ • sin0 m _i ifc )cZ0 m _i, fc 
a m 2 



n+l 



— J l/+ d^ i ( cr fc sin 9i,k ■ ■ ■ sin6> m _ 2 Wa m + ot" m -\ 



(cT fe sin 6>i )fe • • • sin 6> m _ 2 ,fc) 2 ( \/ " m + or m _ x 



In the last step we applied formula (4.11 1 for a = CTkCe m sin ^i^.-- -sin TO _ 2) fc, b = CTkCt m —i sin^i^.- 
• • sin^ m _2.fe and also considered that 

sin(/3 cos a;)(sin a;) iy+1 J v (asmx)dx — 0. 

The integration with respect to the variables 9 m -2,l^ $m-2.n+i follows similarly by applying 
again (4.11 1 and yields 

d0 m -2A • 



r d e m _ %1 . . . r de m -2, n +i TT ^^-^s^-*^^,*™^-,, sin ^+ 

Jo Jo fe=1 



d-m+3 
2 

-2,fe 



— J v+ d^n (cT fc sin 0i, k ■ ■ ■ sin 6 m - 2 ,k\J + a^-i) 
(cr fe sin 6»!, fc • • • sin0 m _ 2ife )i (yo^ + a m _i 



fc=i 



TT \ J v | d-m+i ( cr fc smfi, fc • • smf m _ 3 , fc ^/a^ + + a^_- 



CTfe sin 6x t k ■ - smdm-^k (\[°% 
Continuing in the same way we obtain that 



ra "T" ^"m— 1 ^m— 2 



i^(n,..,r n ,aj = <|2^- 1 r(i, 

Then we can write that 

r((n + l)(2i/ 



^ 1 «+i "+ 1 J, + f _i(cr fc ||q m ||) 
M (cTfcUaJir+f- 1 



(2.8) 



n+l 



(T(2v + d- l)) n+1 i (n+l)(2^+d-l)-l 
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k=l 



dn dr 2 ■ ■ ■ / " dr n TT < +§ " ^- 1VJMI 7 m " y , (2.9) 



In order to work out this n-fold integral, the result (4.8 1 assumes a crucial role. Indeed, we 
apply recursively the formula (4.8) to calculate each integral with respect to the variables Tj. 
Therefore, in the first step we have that 

]Cfc— i r k d d n n 

/ taw* in^+^2 T " +5 ' 1 ( t ~IZ T fe)" + ^ lj ^+f-i( cr "llgmll)^ + |-i(cft-^T A ,0l|Q m 

Jo wmmW) k=1 - - fc=1 

= I ( c| la m ll^+^ (CT " l|a " lir+? " 1(c(t ~g Tfc)ll ^ lir 

n 

^+|-i(cr„||a m ||)J !/+ |_ 1 (c(t- ^r fc )||a m ||) = {cT n \\a m \\ = y) 

k=l 

Mt-T.^ r k )\\a m \\ n-1 
/ ^+.-l (c(t _^ rfe )||« m ||_ yr +^l 

Jo 



(c||a m ||)^-3 



fc=i 

n-l 

J 1/+ |_i(|/)J I/+ |_ 1 (c(t- ^r fc )||a d || -y) 
fc=l 

The second integral is given by 

, 2 



i (r>+w r n: . 

ar„_i 



(c||o; m ||) 6 ^+ 3d — 5 V^fr(2^ + d- 1) 

n— 1 n— 1 

fe=l fc=l 

= (cr„ I |a m 1 1 = y) 

1 ' * dy 



(c\\a m \\)^+^V2nT(2u + d-l)J 

n—2 n— 2 

fe=i fe=i 

(i> + ^)) 3 (cft-E- 2 T fc )||ajr+f^ '- 2 
~(V§^r(3i/+§(d-i)) (c||« m |ir +3d - 4 J 3, +fd - 2 w Z, 7 *^* 

Then, the last integral becomes 

(r(, + ^)f i 



(V^)"-^^ + § (d- 1)) (c||a m ||)(«+ 1 )( 2 ^- 1 )- 2 7 ^ (2 ' 10) 
= (cn||a m || = y) 
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r (v 



d-1 



))' 



<=t||«JI 



{^) n ^Y{n V + ^{d- 1)) (c||a m ||)("+ 1 )( 2 "+ d - 1 )- 1 7 ^ 

(ctlls™!!)^^^" 1 



(r(, + ^))" +1 



(V^) r T((n + l)(^+^)) W|a 



(n+l)(2i/+d-l)-l (ti+1)0+ 



Therefore, plugging the result (2.10) into the expression (2.9 1, and by observing that 



= V^2 



and 



\ = /^)l-(n+l)(2„+d-l) + 1)(2" + ~ 1)) 

y r(^(2i/ + d-i)) 



some manipulations lead to (2.5 1. 



□ 



Let us denote by x m = (xi,X2, ■■■,x m ) & n d dx m = (dxi, dx2, dx m ). Now, by means of 
(2.5), we derive the explicit probability distribution of the random motion (t) , t > 0, which 
is concentrated in the hypersphere H™. 

Theorem 2. The probability law of X m it) , t > 0, is equal to 

P{Xt(t) e dx m } 



II™ 1 

r(^(2i/ + d-i) + |) (c 2 t 2 



-(2v+d-l)- 



r(=±i(2l/ + d- 1) - f + |) 7r m/2( ct )(™+l)(2 i ,+rf-l)- 



(2.11) 



with 



< ct and d > 2. 



Proof. By inverting the characteristic function (2.5 1, we are able to show that the density law of 
the process X TO (i),f > 0, is given by (2.11 1. Indeed, by passing to the hyperspherical coordinates, 
we have that 



?Cd( x m' V 



1 



(27T)' 

1 



dai • • • da ri 



dp / dOi 



de 



m-2 



exp{—ip(x m sin 9% ■ ■ ■ sin 9 m -2 sin </>+•• 

2 4i( 2 ,w-i)4 r( u±i (2 , + d _ i) + i) 



(ct P ) 



i±i(2t/+<i-i)-i 



d^sin" 1 0i • • • sin0 m _ 2 
• + X2 sin 0i cos 02+^1 cos 0i ) 

J^(2v+d-l)-\( ct P)-> 



2 ^ {2 „+d-l)-l n+ l 1 

— — — r ——(2v + d - i) + - 

(2tt)™/2 ^ 2 V ; 2 

2 =±i (2 , y+d _i)_i r ( a ±i( 2l/ + d _ 1 ) + i) 

(2^)™/ 2 ( C i)^(2^-i)-i||x m ||f- 



P (p||x m ||)f-l (cip) H+l(2„ +rf - 1 )-| " 
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r(2±i(2l/ + d-l) + |) 



7r m/a( rf )(n+l)(2«+«I-l)-l F(2±i(2l/ + d - 1) - f + §) 



-(2i/+d-l)-J 



In the first step above we have performed calculations similar to those leading to (2.8) and then 

f-271 



de 1 
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m-2 



II 



> sin™ 2 6»i • • • sin 6> m _ 2 



(2.12) 



exp {—ip(x m sin #i • • • sin Q m -i sin < 
J™ xMlxJI) 



£2 sin 0i cos 02 + a^i cos^i)} 



(27T)" 



(pIKII) 1 



while in the last step we have used formula (4.101 for v = ^i^-(2u + d— 1) — |, fx ■■ 
and b = ||x„, II. 



1, a = ct 
□ 



Remark 2.1. By taking into account (2.11 



observe that X^ (t) , t > 0, represents 



isotropic random walk; that is its probability distribution depends on the distance, from the ori- 
gin, of the position of the random walker. Furthermore, by setting v — in (2.111, we reobtain 
the result (2.26) in De Gregorio and Orsingher (2011). 

Remark 2.2. The cumulative distribution function for the process X m (t) , t > 0, is equal to 



~ m 

p{*L(t) g n™} = / p^ d (* m , t) n d Xl 
j ^ t=i 



2T(a± 1 (2i/ + d-l) 



r(2±i(2i/ + d-l) 



l)r( ! 



„ m _x(c 2 * 2 -p 2 ) 



l -(2i/+d-l)-- 







^(f»+l)(2v+d-l)-l 



-dp 



(2.13) 

" ■- J ■••->" - - 1 - ■ " = g € N, we can write down P{X^(t) G H™} in an 
alternative form involving a finite sum. Then, one has that 



withQ<r<ct.If r ^{2v + d-l)- ., 



2r(a 



1) 



2^ 



r(q + i)r(f) J Q 

T(q + f + 1 



p 

(r/ct) 2 



dp 



r( o + i)r(f)y 



-! (C^-P 2 ) 

(ct) 2 i+ 



"D-D 



2 ; fe=0 



y Vet 



2/H 



Remark 2.3. 5y assuming that the number of steps n + 1 is a random value, it is possible to 
get the unconditional distribution for X m (t),t > 0. In other words, we suppose that the number 
n of changes of direction are governed by a fractional Poisson process N%(t),t > 0, introduced 
by Orsingher and Beghin (2009). Therefore, we have that 



P{NX(t) = n} 



(At)" 



r(=±i(2!/ + d - 1) + i) E u+ *-i „ +4 (At) ' 



n = 0,1,2, 



(2.14) 



, x G K, a, /3 > is the Mittag-Leffier function. Furthermore 



where E a ^(x) = J2k=o k\r(lk+p) » • 
N%(t),t > 0, is supposed independent from r and 0. Then, in order to obtain the unconditional 
distribution o/X m (t), £ > 0, we can average (2.111 wii/i i/ie distribution (2.14). 
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Remark 2.4. Let us consider the radial process R^it) = ||X_f n (t)|| = Vz^Li Xi{t),t > 0. We 
observe that 

P{R d m (t) <r} = P{X? n (t) E K n }, 



withO < r < ct. Therefore, letf£ m d( r i~t) ^ e ^ e density function ofRf n (t), the result (2. 13 \ implies 
that 



fn,m,d( r ' 



2r(=±i(2l/ + d- 1) + i) r m-l( c 2 t 2_ r 2)2±l(2„+d-l)-^ 
r(=±i(2l/ + d- 1) - I + §)r(f ) ( c< )(n+l)(2„+d-l)-l 



Furthermore, by taking into account (2.15), we derive the moments o/i?^(i), that is 

E « (t)F - im* + i-n + i)im 

with p > 1 . 



r(2±i( 2l/ + d_i) + E±i)r(s) 



(ct) 1 



(2.15) 



3 Random flights with v = 1 

As observed in the previous Section for a general value of v, we are not able to work out the 
integral ( |2.1[ ). Nevertheless, it is possible to derive a closed-form expression for the characteristic 
function of X rf (t),i > 0, for some values of v. The next Theorem provides the explicit char- 
acteristic function for the random flight with v = 1. It is worth to mention that this result is 
remarkable since in this context it is hard to obtain explicit results. 



Theorem 3. For v = 1, the process X d (i), £ > 0, admits the following characteristic function 

n+l-j 

/ i \ I " d >."t-l; 



/7rT((n + l)(d + l)) 



. (d+l)(n+l)-l 



J=0 



n+ 1\ \\\a d \\ 2 2 



J („+i)(d+3)-(2j+i) (ct\ |a d | I) 



r( 



(rt+l)(d+3) 



i) (c*||a d | 



(,i + l)(d-l) + 2j-l 



(3.1) 



imtt n > 1 and d > 2. 



Proof. In order to prove the result (3.1l, we will use the same approach and tools exploited in 
the proof of Theorem [l] Nevertheless the proof of (3.1 1 is not a simple adjustment of the proof 
of Theorem [T] and it requires a particular care as we will show in the next steps. 

We start by taking into account the expression (2.1 1 for v = 1. Thus, we have to handle the 
following integral 



X n (n, ...,r„,a d ) 



r(l + d/2) 



n+l 



Tt + 1 



dO, 



d-2,1 ■ 



d6 d -2,n+l / #1 ' 
JO 



(3.2) 

2n 



k=l 



JJ "I exp{icrfc(o!dSm6i,fc sin0 2 ,fc • ■ • sin0d-2,fc sin^ + sinfli,* sin0 2 ,j • • • sind d -2,k cos0 fc 
a 2 sin 6» lifc cos 2 ,fc + a x cos 0i,k)} sin 9% k 



sin 3 6» d _ 2 ,fe sin 2 fe 
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For the n + 1 integrals with respect to 4> k , we get that 



2tt j-2n n+1 

#! • • • [ 



(i 



+ 1 e ic1 "fc sin 01,* sin 9 2 ,fc---sine d _ 2 ,fe(Qd sin fc +Q d _i cos fc ) ;2 



sin <p k a<p k 



n+l ( r 2TT 

n 



*:=i 



g icr fc sin 6li ifc sin 2> * "•sin 6 d -2,k (ad sin ^fc+aj-i cos fc ) g j n ^2 



"+ 1 f Ji(cr fc sin^i,fc sin 2 , fe • ■ • sin0 d _2,/cA/a d + «d-i) 

= (2vr)" +1 J] 7 

fc=i I CTfe sin 0i,fe sin 2 ,fe ••• sin 0d-2,fe a/ a d + 



+ a d-i 



J 2 (cr fe sin 6»i !fc sin 6» 2 , fc • • • sin0 d _ 2 Wa d + a^) >, 



where in the last step we have used the result (4.3). Therefore the integral (3.2) becomes 
^(n, -,T„,a d ) 

r(i+d/2) \ + (27r)n+1 r d0ii ... r dein+1 ... r d9 d _ 2S ... r d6 d _ 2 , l+1 

7T2 / Jo JO JO JO 

n+l r 

J~J J e *cr fc (a d _2 sinej.fc sin 6 2 ,k •••cos 6» d _ 2 ,fcH ha 2 sin0 1:fc cos0 2 ,fc+a:i cos6»i_ fc ) g j n Qd ^ _ _ g j n 3 ^ 

fc=l I 



Ji(cTfcsin0i !fc sin0 2 ,fc • ■ • sin0 d _ 2 ,fcWa d + a2_ x ) 



cr fc sin 0i ife sm 2 ,fc • • ■sinyd-2,fc\/a d + a d _ 1 



or, 



-J 2 (cT fc sin 0i >fe sin0 2 , fc • • • sin0 d _ 2 Wa d + a 2 d _ 1 

"d -r 

The integrals with respect to 0d-2,k are equal to 



pit pn n+i [ 

/ ^d-2,1 • • • / dfl«j-a,n+i II i 

- 70 - 70 fc=l I 



icT k (a d - 2 sin 0i >fc sin 2 ,fc • -cos 9 d -2,k) a 



Ji (cr fe sin 0i jfe sin 2jfc • • • sin 0d-2,fc \/ + «d-i) 



cr k sin0i, fc sin»2,fc • • • smt^Wo^ + a d _ 1 



n+l ( pTr/2 

n 2/ ^ 

s:=i i j ° 



-J 2 (cT fc sin 0i jfc sin 2 , fc • • • sm9 d _ 2 ,k\ u d + «d-i) 



-2,fe cos(cr fc ad_ 2 sin 6 x<k sin 2 ,fc • • • cos 9 d -2,k) sin3 #d-2,fc 



Ji(cT fc sin0i ifc sin0 2 ,fc • • • sin d - 2, k \ atj + a^-i) 



CTfc sin0i jfc sin0 2, fc • • • smtfd-2,fc\/Q!d + "d-i 



-J 2 {cT k sin0i ifc sin0 2 ,fc • • • sin 9 d ~2,k\ a d + oq^ 



l d "d-i 
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n+l 



J3/2 (cTfe sin 9 1 tk sin6» 2 ,fc • • ■ sia9 d -3,k\ / oq + oq_ x 



l d-2l 



n 2 

k=i { V Z L (cr fc sin 6»i^ fe sin B 2 . k ■••sin 6» d _ 3 , fe Jc^ + a 2 ,^ + c^_ 2 ) 3 / 2 
"d' / 5/2(cT fc sin6'i ifc sin6'2,fe • • • <smQd-z,k\j a \ + a d-i + a d-2) 
(cr fc sme^k sin 6 2 ,k • • • sin0 d -3,fc) 1/2 ( / a 2 + Q ,2_ i + a 2_ 2 )5/2 



where in the last step we have used formula (4.11 ) for 



a = CT k smU ltk smU 2 ,k • ' • smf d _ 3 , fc y a 2 , + a 2 , _ x 

6 = cT fc a (; _2 sin sin 6* 2 ,fc • • • sin 6» d _ 3 , fc 
and ^ = 1,2 and also considered that 



(sinx) sm(bcosx)J u (asiiix)dx = 



By means of the same arguments, we can calculate the further integrals with respect to 
0i,fc) •••) l,fe 5 j @d-3,k and then we obtain that 

n+l 



l 1 n {r u ...,T n ,a d ) = {2 d / 2 r(l + d/2)} J[ 



™ +1 TT / J d/2( cr fcll«dll) a 2 d J d/2+1 (cT k \\a d \ 



ifil (cr fe ||a d ||)^ ( CTfe)rf /2-i| K || d /2 + i 
The characteristic function for ^ = 1 becomes 



T (3.3) 



J 7 . 



l(a d ) = {2 d / 2 r(l + d/2)} 



Tt + l 



t-Tl 



(Zto • • ■ 



<^7"n — jT n ) 



n+l 

n 



/ ^rf/2 (en- 1 |a d 1 1 ) a 2 Jd/2+1 (cr fe I |a d I 



i I {cr k \\* d \\) d/2 (cr k ) d /^\\a d \\^ 



(3.4) 



By developing the product appearing in (3.4 1, the characteristic function Fni&d) reduces to 
a sum, over to 2™ +1 elements, with generic terms given by 



C 3 B, fe) 



(3.5) 



= f 2 d / 2 T(l + d/2) \ n r((n + l)(d+l)) (_ al 
T{d+l) J t {n+i)(d+i)-i y \\a d 



2 \ n+l- j 
2 



rfdn / r 2 d dr2 • • 



t-Ti 

ll Jo 

o 



T ni-l^ T ni-l 



n 



Jd/2(cT k \\a d \\) 



\ (cr fc ||a rf ||)^ 



T rii + l^ T ni+l 
T n 2 -1 UT ™2-1 



r ni+2^ r ni+2 ' ' ' 



Tfc ■ /-'-E^ , J* J d / 2+ i(cr fc ||a d ||) 



a« 2 n 



fc— ni + l 



(cnbllfidH)"/ 



2-1 



*-E fc =i T fc 



T n s - 1 +2^ T n 3 -l+2 ' 



14 



Jo ^-1^-1 ^ r ns dr n3 (ar k \\a d \\)^ 

/ r n a + l dT n s + l / T n s +2 d Tn s +2 ' ' ' 

Jo Jo 

with < ni < «2 < ... < n s -i < n s < n + 1, j = ni + JI3 — ri2 + ... + n s — n s _i and < j < n + 1. 
Therefore, we focus our attention on the calculation of (3.5 1. 
We deal with the following (n — n s )-fold integral 

=/ T n 3 + ldTn s + l I T n B +2 dT n s +2- ■ ■ 







Jo Jo ,„_, ,„ 



Jd/2+i(cT k \\a d \\) 
, „. +1 M^ll)^" 1 



by applying formula (4.8 1. Indeed, the integral with respect to t„ is given by 
t-E^-r* , " . J d / 2+ i(cT»||a d ||) Jd/2+i(c(t-ELi r fc)Hgrfll) , 

T ife (cT-nll^ll)"/ 2 - 1 Wt-ELi^)!!^!!)^ 2 - 1 T " 

rt-E^i^ (cr n ||a d ||) d / 2 + 1 (c(t-ELir fe )||a d ||) d / 2+1 



J o ( c || Qd ||)2d Jd/2+i{cT n \\a d \\)J d/2+1 {c(t- 2_^T k )\\a d \\)dT n 

(cr n 1 1 oj] I = to) 

Mt-T^l^WsuW w d /*+\c{t-Y!-{ r k )\\ a d \\-w) d / 2 + l 7 ,. T ^ ... 



1 r 2 ((d + 3)/2) 



fc=i 

n— 1 n— 1 

(c(t - E ^)ii« d ii) d+5/2 ^ + 5/ 2 ( C (i - E ^)ii« d ii) 



WKII) 2d+1 ^r(d + 3) 



where in the last step we have used formula (4.8 1. The integral with respect to r n _i becomes 

1 r 2 ((d + 3)/2) 

(c||a d ||) 3d+1 V2^T(d + 3) 

/ \cr n ^\\a d \\) d / 2+ \c(t-Y,r k )\\^ 

Jo k=i k=i 

= (cr„_i||a d || = w) 

= 1 r 2 ((rf + 3)/2) 
~ (c|KII) 3d + 2 V2^r(d + 3) 

r c (t~ E n — 2 r fe) Q II n— 2 n— 2 

/ ' * ^ d / 2+1 (c(i - ^ ^)||« d || ~ ^) d+5/2 ^/ 2+ i(u;)^ +5/2 (c(* - ^ ^Oll^ll - 

Jo ._, 



k=l k=l 

+ , ,. ^ 2rf+4T , ^ ... ... 

T-fc)||adll), 



, ,, ^rr, T ((d + 3)/2) _ y 1 d+4 _ y 

(c\\a d \\) 3d + 2 (V2^) 2 T(Ud + 3)y V ^ll-dll^ 5 <i+4V V 



fc=i fe=i 
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where in the last step we have used the formula (4.8) again. We can continue at the same way 
and then we have that 

r^r \ = 1 r"~"-((d + 3)/2) 

(c||a £i ||)( d + 1 )("- n =+ 1 )- 2 (v^)T,-n.-ir(2^U(d + 3)) 

ft— Y] n t Tfc n a +l 

/ = dT. ns+1 {cT ns+1 \\a d \\) d l 2+ \c(t~ J2 ^)IKII) < "^ 2,(d+3) " 1 

J ° k=l 

n s + l 

Jd/2+l{CT nr +i\\c^\\)J ( n - n M4+3) i (c(t - T fc)ll«dll) 

fc=l 

= (cT„ 3+ i||a rf || = w) 

1 r™-"=((d + 3)/2) 



(c||a d ||)( d + 1 )("-™^+ 1 )- 1 (v^)n-Tu-ir(5^(d + 3)) 

("-"s)( ' i ■ ■ i 

r fc)||a d || -^J 

k=l 

n s 

Jd/2+lMJ (n-„ s )(d+3) i (c(t - V, Tk)\ \Qd\\ ~ w ) 

2 2 £ ' 

fc=l 



2 

/ dW/^i - Tk)\\*A \-wy 

JO 



= (c||a d ||)(^)("-"-+i)-i(v^F)"-"-r(n±kzn.(d + 3 )) J ^^ -| ^ ~ Ib**)^ 
Now, we perform the calculations concerning the (n s — n s _i)-fold integral 



l_! + l dT n s -l + l / r n s „ 1 +2 dr " s -l+2 ' ' ' 

Jo 

By means of the same approach adopted above, we can write down that 

r" +1 -™=((d + 3)/2) 



(c||a d ||)( d + 1 )(«-«=+ 1 )- 1 + d (V2^)"-^r(^±^(d + 3)) 

rt-J2k!=i 1 T k ™ s 

/ ' dr„ s ( C r n JK||) d / 2 ( C (t-^r fe )| 

- 70 fe=i 

Jd/2(cT n , n«^ii) ■ /(»+i-n.)(j+3) i (c(f - Tfc 



2 

•£^fc=l Tfc , , , , (n+l-n s )(d+3) 

:::<yi) 
fe=i 

(cT„J|a d || = w) 

r n+1 -"-((d + 3)/2) 

(c||a rf ||) (d+1)( " _ " s+1)+d (v^i : )"- n »r(5±^(d + 3)) 



dtim; d / 2 (c(t- J2 T k)\hd\\ -w)~ 



k=l 



1G 



n 3 -l 

Jd/2(w)J (n+i-n e )(d+3) i(c(t- ^2 r k )\\a d \\ - w) 

k=l 

r" +1 -"s((rf + 3)/2)r((rf+l)/2) 

( c ||o; d ||)(rf+i)("-n 3 +i)+d(y27r)«+i-K ar (- ("+ 1 -(»»- 1 ))( d + 3 ) _ ^ 

(„ + l-(„ a -l))(d + 3) 3 
I 2 2 < 

fe=l fc=l 



( c (*~ 51 T k )\\a d \\) in+1 '"^ 1)><d+3> 2 J („+i-(„ 3 -i))(d +3 ) 3 (eft- T k)\hd\\)- 



Then, carrying on the same calculations for the integrals with respect to the variables 
Tn a -i, -,t„ s _ 1+2 , we get that 

r n+i-n 3{{d + 3 )/2)r».-".-i- \(d + l)/2) 



(c||g d ||) (d+1)( "~"- l) ~ 1+d (\/2^)"-"- 1 - 1 r( ( "~^- l) (d + 3) - (n 5 -n s _! - 1)) 

/•t-Z;"! -1 Tfc n s _i + l ( _ ^ 

/ " dr„ s _ 1+1 (cr„ 8 _ 1+1 ||a d ||) d / 2 (c(t- £ r fe )||a d ||)^ Z ^( d+3 )-^-"- 1 -^ 

n e -!+l 

J d /2(cr ras _ 1+1 iiQ d i|)J ( n-n 1} _ (eft- ^ r fc )||a d ||) 

fc=l 

= (cT„ s _ 1+ i||a d || = w) 

r"+ 1 -" s ((d + 3)/2)r n ' s -™ s - 1 - 1 ((rf+l)/2) 
~ (c|| gd ||)( d + 1 )("-^-i+ 1 )- 1 (\/2^)"~" s - 1 ~ 1 r( (n '^- l) (^ + 3) - K - n s _! - 1)) 



/ ^/ 2 (cft-£r fe )||a d ||-™ 



( " " 2 3 - l) (d+3)-(« s -« 3 _ 1 -^) 



J rf /2HJ ( n-, 1} wa _ ws _ i , (Cft- £ -«o 



fc=l 



((rf + i)/2)«+ 1 -« 3 r™+ 1 -"- 1 ((d + l)/2) 



~ (cila d |l)("+ 1 -^-0( rf + 1 )- 1 (\/2^)"-^-ir( ( " +1 "" s 2 - l)(d+3) - (n s - n._i)) 

( n+ l_ n )( d+3 ) 

(c(t- £ r fe )||a d ||) » ("--"-i +1 / 2 >J ( „ +1 -„ l)(d+ 3) (w< ras _ l+1/2) (cft-^ ^fel 



fe=i fc=i 

where in the last step we have also exploited the formula T(l + z) = zT(z). 

Now, bearing in mind the same steps used to work out IC ns (a d ), we are able to explicit, 
omitting the calculations, the following integral 



= / ^3-2 + 1^-2 + 1 / T n s - 2 +2 dT n s - 2 +2 ■ ■ ■ 

Jo Jo 

r^^ T \ d dr , r E ""'\^ dr K7-»»>(a) TT J W^IKII) 
Jo T '-- 1< - 1 - 1 io ^-A—K ^Jii^ (cTfc ||a d ||)^-i 
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((d + i)/2)"+ 1 -« 3 +«-i-«-2r™+ 1 -"»- 2 (((i+ l)/2) 



(c||a <i ||)("+ 1 -"-=')(' i + 1 )- 1 (v^i : ) n -"- a r( ( " +1 ~ n ' 2 - 2)( ' i+3) - (n 8 - n a _i)) 

n,_2 "s-2 
(C(t- E^)il^ll) (n + l "" S2 " 2)(<i+3) " (raa "" S - 1+1/2) J (n +1 -n s _ 2)(d+ 3) K n ._ i|1/2) (c(*- E^H^ID- 



k=l ~ fe=l 

Analogously to /C n * _1 ' ns (a d ), we get that 

£n.-s,n._ 2 ,n._i,n.^ 



= T n s - 3 + l dT n s - 3 + l T n 3 _ 3 +2 dr n 3 - 3 +2 ' ' ■ 

Jo Jo 



/ <_ 2 -idr„._ a _i / <_ a dr n ._ a X7 l —"—"-(a <I ) JJ 

• A ' J ° " fe=n s _ : 

((d + 1) /2)"+ 1 -™-+™»-i-™-2r™+ 1 - ns - 3 ((d + 1) /2) 
~ (c||a d ||)("+ 1 -"--3)^+i)-Hv^) n -"-- 3 r(^±^ -- 3)( ' i+3) - [K - + K-2 - n a - 3 )]) 

n s -3 

(n + l-n s _ 3 )(<i+3) ,,,,„, 

(C(t- E r fe)l|adll) ~ 2 [(n.-".-i) + (n.-2-n.-s)+l/2] 

fe=l 

x j ( „ +1 -„._ 3)(d+ 3) _ [(w ,_ n ._ l)+(w ,_ 2 _ w ,_, )+1/2] - E ^H^ID- 

Finally, by applying recursively the same arguments exploited so far, we get that 

-/;r ns y 

f 2 ^r(i + d/2) )" +1 r((n + i)(d + i)) / o? d y +1 ~ j 

\ r(d+l) J t(n+l)(d+l)-l ^ Hodiiv 

Jo Tl ^ Jo ^ " ' io J ^dr ni [[ {cTklhdll)d/2 

/ T ni+1 d,T ni+1 T ni+2 dT ni+2 ■ ■ ■ 

Jo Jo 



"s-2 



Jd/2(cT k \\a d \\) 



/ <-i^n 2 -i / <dr„ 2 ] 

Jo Jo , t 



drf/2+i(cT fc ||a d ||) 
, _„, +1 (c^HaalD^-i 



((d + 1) /2)™+ 1 ~[( ns_ ™ s - 1 )+(™ s - 2_ ™ s - 3 )+---+(" 3_ ™ 2 )]r" +1_ " 2 ((d + 1) /2) 



( C || gd ||)("+i-" 2 )(^+i)-i(^)"-" 2 r( ( " +1 - r 2 2)(d+3) - [(n s - n s _!) + (n s _ 2 - n s _ 3 ) + ... + (n 3 - n 2 )]) 

™2 

(C(t -E T fc)llgrfH) ( " + 1 '" 22)(d+3> - [(n.-n.- 1 ) + (n.- a -n._ 3 ) + - + (»3-n 2 ) + l/2] 



fc=l 

"2 



^ ("+i- ? )^+») _ [(w ._ n ._ l)+(w ._ 3 _ w ._3 )+ ... +(n ,_^ )+1/2] (c(*-E 7ifc )ii^i 

fc=l 

f 2^r(l + d/2) )" +1 r((n + l)(d+l)) / ag 

i r(d + i) j t (n+i)(*fi)-i ^ HojHV 
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((d + i)/2)" +1 ^r" +1 ((rf + i)/2) 



(c| \a d \ |)(n+i)(d+i)-i (v^)»r( (" +1 H rf + 3 ) - j 

V5rT((n+l)(d + l)) (-Hpt^ 1 )/ 2 )) 



... (n + l)(d+3) ■ 

(cfc||) ' J ~ 



i±t21_y + i)(c*ll«dll) 



n+l-j 



2 (d+1)( ; +1) ' 1 r( 



(n+l)(d+3) 
2 



(cfcll) 



(Tl+l)(d-l)+2j-l 



- J (^.+i)(d+3)-(2j+i) (ct||q rf ||). 



The integrals J^ 1 ^ 2 ' ,ns {OLd) have different configurations involving the Bessel functions 
Jd (crfc || a d \ |) and Jd_ 1 (cTfe||a d ||), according to the values m,n2,.... 



previous calculations, we observed that the value of J 7 



n\ ,ri2 ,. . . ,n £ 



Nevertheless, from the 
(a d ) is the same if the j does 
not change for different combinations of Ui,ri2, ...,n s . In other words, in order to calculate the 
values of the different terms of the sum appearing in (|3.4|), we have to consider the number of 



Bessel functions with index = (or f — 1) appearing in (3.5 1. Then, we can conclude that 



n+l 
3=0 



n + l 
3 



7" 



and the proof is completed. 



(3.6) 



□ 



Remark 3.1. The analytic j or m of the characteristic junction T\{ol^ is more complicated than 
its counterpart in the uniform case (y = 0) (see formula (2.1) in De Gregorio and Orsingher, 
2011) which is only given by a Bessel function (with a suitable constant). This is because the 
random motions considered in this paper have drift. Indeed, the choice of directions non-uniformly 
distributed, for each step performed by the motion, implies the loss of spherical symmetry. 



Remark 3.2. By means of (4.4 1 and (4.5 1, the characteristic function of~K d (t),t > 0, could 
be calculated for v = 2,3. Clearly, in these cases the necessary calculations for completing the 
proof follow the same steps of the proof of Theorem [3| are very long and cumbersome. Indeed, 
the expression (|3.4|) involves the product of linear combinations of Bessel functions. 



Theorem [3] allows to show our next result. We are able to provide the density law of X d (t), t > 
0, (up to some positive constants). 

Theorem 4. The density function ofX. d (t), t > 0, for v = 1 assumes the following form 



n+l 



r((n + l)(d+l))_y, f _ 1)n+H 



7T^(2rf)(«+ 1 )( d + 1 )- 1 ^ 
k 



h rf^-fc + i 



3- 

2k 



n+l 
3 



r( 



(w+l)(rf+3) 
2 



j) 



(c^-lfel 



'—k 



(3.7) 



where ||x d || < ct and ak, n +i-jS are positive constants defined as in (4.7l. 



Proof. We invert the characteristic function (3.1l as follows 



(2 



Jr* f = \ 



/TrT((n + l)(d+l)) 2±i 



(27r) d 2 



(d+l)(r» + l)-l 



3=0 



n + l 
j 



(n + l)(d-l)+2j-l 

(ct) 2 r( 



(n+l)(d+3) 
2 



■i) 
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2 \ n+l-j J (r 



+ l)(d+3)-(2j + l) 



(ct\h d \\)<L 



/7rT((n + l)(d+l)) 



n+1 



(27r) d 2 



(d+l)(n + l) 



_^£ ( _ ir+1 - 



'n + 1 



(n+l)(d-l) + 2j- 



j — da. 



3=0 



- !— ! - Lj -.' 1 ^/ (n+l)(d+3) _ ^ 



J / (rf) 2 J r(* 



(l Jo ^0 ^0 



2 71 



2 / ^^g — 'ip(^d sin f?i ••• sin 0+a; c i- 1 sin ■•• cos 0+. . .+X2 sin cos ^2+^1 cos 0i) 



(sin 0i • • • sin Od-2 sin < 



„, . J (n+l)(d+3)-(2j + l) (dp) 

x2(n+l-j) 2 „;„ a d-2 



(n + l)(d-l) + 2j-l 



sin i 



■ sin 0d-2- 



By taking into account the formula (4.7 1 , the integral with respect to <j) becomes 

e -ip(x d sin di — sinfl d _ 2 sin0+a: ti _ 1 sin0i — sine d _ 2 cos 0) g m 2 (™+ 1 -i)^0 



n 



n+1- 



2. J; 



j (-^) k ak, n +i-jxf J n+ i- ]+ k(psiiie 1 ■ ■■sm9 d _2)Jx' 2 d + x 2 d _ 1 ) 



fe= o 2 fc (psin6>i • ■ ■sin0 d _ 2 ) n + 1 -i- k (^Jx 2 d + x 2 d ^ 1 ) n + 1 -^+ k 



Then 



Pn(x<i»*) 



n+1 



y^r((n + i)(d + i)) g^^n+i-j - 1 



(27r)d -l 2 ^±^ 



J=0 



'd+l^+i-j 



j /( c t) (n+1)( V )+2i - 1 r( (" +1 » d+3 ) -i) 



" + 1 i ( -1 \fc„ „2fe /-co 

2 k 



J (Ti+l)(d+3)-(2j + l) (dp) 
1 



(n + l)(d-l) + 2j- 



E 

fe=0 

/ d6\ ■ ■ ■ I d6 d 2 g^V( a; ti-2 sinSi- - sine d _ 3 cos6» d _2+ - +2:i cosSi) g - n ^d 

Jo Jo 

( S infl 1 ---sin0 (i - 2 )" +1 - j+ fc 

(V^ + *Li) n+1 - j+fe 



d-2 



• • • sin W d _ 2 



-J„+i_j+fe(psin^i ■ ■ ■ sm9 d -2\ x 2 d + x 2 ^) 



By means of the same approach used to get (2.8 1, the (d — 2)-fold integral with respect to 
the angles becomes 



/>7T nTT 

j dOl ■ ■ ■ / dOd 2e~ ip ( Xd ~ 2 sin ^!'" sinfld-3 cos 8 d - 2 + ---+xi cosfli) g j n gd -2 

Jo Jo 



■ sm 6» d _ 2 



n+l-j+k 



(sin 0i • • • sin0 d _ 2 ) 
(2 7 r)f- 1 



-J n+ i_j +fe (psin^i • • •sin0 d _ 2 Wa;3 + 



pf- 1 ||x d ||"+ 1 --''+ i: + 



iZT J n+i-i+fc+|-iWi^dlD- 



Therefore, the formula (4.101 implies that 
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Tl+l 



0fr((rc+l)(d+l))^ ( 1)n+1 . 



, , d (d+l)(r,+ l)- 

(27r) 2 2 



3=0 



71+1- j 



(— l) fe afc,)i+i-j 



2' 

(by (pOt) 



,.2/r 



i 1\ /rf+l\«+ 1 -J 

OO J („ + l) (d+3 )-(2j + l) (ctp)J Ti + 1 _ J - + fc+ ^_ 1 (p||x w ||) 



|7l+l-j + fe+|-l 



n+1 



0fT((rc + l)(d+l)) v , ,„ +1 - 
(2tt) 2 2 2 a . =0 

I, - J-J "/,M: J- ; 

fc=0 



n + 1 



d + l\n+l-j 



3 /(ct) 

J2h 



(n + l)(d-l)+2j-l 



r ^ (n+l)(d+3) _ 



2^ 



n(d+l) 
2 



fc + 1 (ct) 



(n + l)(d + 3)-(2j + l) 



(a 2 -||x d || 2 )^ 



n+1 



r((n + l)(d + l)) +1 . 
7r^(2rf)("+ 1 )( rf + 1 )- 1 ^ 



n + 1 



j yr( 



d+1 xn+l-j 



>(ra+lKj+3) 



■i) 



71+1- j 



E 



(— 1) afe ; „-|_i_jX d ' 



2// 



(d+l) 



fc+1 



(A 2 -||x d || 2 )- 



'—It 



dp 



□ 



It is worth to point out again that the result stated in the above Theorem does not provide 
the whole analytic form of the density law of the random flight with v = 1. Indeed the constants 
appearing in the expression (3.7 1 are not determined for an arbitrary value of n, but only for 
some values of n and n + l—j. For instance ao, n +i-3 = (2(n + 1 — j) — l)(2(n+ 1 — j) — 3) • • • 3- 1 
and a n +±-j >n +±-j — 2 n+1 ~i (see formula (4.7 1 in Appendix). Furthermore, in a few cases we can 
explicit (x d ,i). By taking into account (4.6 1, for n = 1 and n = 2, after some calculations, we 
obtain that 



r(2(d + i)) 



tt^ (2ct) 2d + 1 (d + 2)r(d + l)r(^fi) 



( C 2 i 2 -||x d || 2 )^- 2a; 2 ( C 2 i 2 



d-i 



(d+l)x d (c 2 t 2 



and 



pi fx t) = r(3d + 3)(d + i) 

n^(2ct) 3d + 2 T(d-l)T(-^-3)(3d+7)(3d + 5) 
j 4(d + 4) (c 2 f2 || Xd ||2)d+i | 2x 2 (6d 2 + 6d + 8) ( ^ 2 



(d + l)d(d-l) 



(d + l)d(d-l) 



(3.8) 



s^^-ifeHY^ + ^d+i^cA 2 



|2\d-2 



(3.9) 



Let f% d( r 't) ^ e the density function of i?d(0 = ||X rf (f)||, i > 0, with < r < ct. As done for 
obtaining ( |2. 15 > , from (3.8 1 and ( |3.9[ ) we get that 



r(2(d+l))2v^F 



(2ct) M +i(d + 2)r(d + i)r(<*=i)r(f ) 
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d=2 

d=3 
.... d=4 



Figure 3: The behavior of the density functions f\ = fl d {r, 1) (on the left) and fa = fnjir, 1) 
(on the right) for d — 2, 3, 4 and c = 1. 



d-1 



r d - l {c z ? -r 2 )- 



2 r d+1 (c 2 t 2 - r 2 )^ + 3{d+1) r d+3 (c 2 t 2 - r 2 f 



d(d + 2) 



(3.10) 



and 



fL(r,t) = 



r(M + 3)(d+l)2y/ir 
{2cty d + 2 Y{d - l)r(?!f - 3)r(f)(3d + 7)(3d + 5) 

' 4 (^+ 4 ) r d-i fr 2 f 2 _ 2, d+ i , 2(6d 2 +6d + 8) j+lf 2 2 _ 2 d 
(d+l)d(d-l) 1 ' (rf+l)d 2 (d-l) 1 j 
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d(d + 2) 



„d+3^2 t 2 _ r 2^d-l + 40(d+l) _ r 2-jrf-2 



d(d + 2)(d + 4) 



(3.11) 



Remark 3.3. TTie expressions (3.8 1 and (|3.9|) represent proper density functions. Indeed, it 



is possible to verify after some calculations that J^ d p„(x d ,t)dx d = l,n = 1,2. Furthermore, 

p\(x d ,t) and p\(x. d ,t), or equivalently f\ d {r,t) and f2d( r ^)> are non negative functions (see 
Figure^. 

Remark 3.4. It is useful to compare the density functions p\(x d ,t) and p\{~x. d ,t) with the prob- 
ability distributions obtained in the uniform case (v — 0) by Orsingher and De Gregorio (2007). 
As emerges from Table the uniform random flight and the random motion with drift lead to 
quite different distributions. 



Remark 3.5. From the results (3.7l ; (3.8 I and (3.9 I emerge that the probability distributions of 
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d = 3 
d = 4 



2/ = 


1/ = 1 


1 


15 


2-7rct^/c 2 t 2 -||x 2 || 2 


2 4 7r(ct)5 


2+2,11 n 


21 


T(2ct) 2 ||x 3 | 


2 4 ir(ct)7 


2 


3 2 .7.5 


jr 2 (ct) 4 


257r 2 (ct)9 



^ 1 1 (A 2 - iix 2 n 2 )i - xWc^-\\x 2 \\> + 



d = 2 



d = ■ 



{|(A 2 - ||x 3 H 2 ) 2 - 2x|(c 2 t 2 - ||x 3 || 2 ) +4x|} 

(c 2 t 2 - ||x 4 || 2 )S - 2xt(cH 2 - ||x 4 || 2 )i +5xt^cH* - Uxjp} 
Table 1: Density function p\ (x^, t) for iv = 0, 1, and d = 2, 3, 4. 
i/ = 1/ = 1 



1 

T (ct) 2 



2_jr 

T(ct) s 13-ll 



(c 2 t 2 - ||x 2 H 2 ) 3 + ^xl(c 2 t 2 - |x 2 || 2 ) 2 



-2x\{c 2 t 2 - ||x 2 



-2^ 



^(c 2 * 2 - iix.ii 2 ) JlS^ \ - ii^iip)- + -<>^ - 1is.11 3 >' 

-4(c 2 t 2 -||x 4 || 2 ) a + |4(^^-||x 4 | 
Table 2: Density function p%(it d , t) for 1/ = 0, 1, and d = 2,4. 



the random flight with v — 1 are a linear combination of n + 2 functions 

xf(ch 2 -\\^\\ 2 )^- k 



with k = 0, 1, n. + 1. Furthermore, (3.7) permits to claim that, as expected, the random flight 



with drift is not isotropic. It is interesting to point out that the observer in W n ,m < d, perceives 



the random motion as an isotropic walk (see Remark 2.1), while the original motion is non- 
isotropic. 

4 Appendix 

In the random flight problem the integrals involving Bessel functions have a crucial role. For 
this reason in this Section we summarize some important results which are often used in the 
manuscript . 

In particular we focus our attention on the following integral 



Jz(a cos 9+/3sine) q^q 



(4.1) 



with a,/?£l and v > 0. We are able to calculate (4.1 ) in some particular cases. For example, 
for v = it is well-known that 



el z( a c OS e+p S in0) de = 2nJ (zyJa 2 + p 2 ). 



(4.2) 
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If v = 1, we get that 



f * ^(acos^sin*) sin 2 0d9 = 2w Jl(z Va 2 +P 2 ) _ _P_ {z ^^2 } I _ (43) 

as proved in De Gregorio (2010). The previous result has been obtained by expanding the 
exponential function inside the integral. Hence, one observes that splitting the integral into two 
parts, i.e. J Q 27r = J* + jj^, after a change of variable, the sum is not equal to zero if the index 
values are even. 

Following the same approach developed in De Gregorio (2010) we can write that 

/■2-7T 

/ exp{iz(a cos 9 + f3 sin 9)} sin 4 9d9 
Jo 

= 4 y ^tz 2k y ( f) a 2r (3 2 ^ f 72 cos 2 - 9 sin 2 ^"^ 2 ) 9d9 
^ o (2fc!) 2^\2r) J 

_ 9 f - til. J»* V ( 2k \ ^ fl2(*-r) r(r + l/2)r(fc-r + 2 + l/2) 

- 2 ^ n (2*1) hX^) 1 r(fc + 3 ) 



k=0 v y r=0 

oo , , .1. ot. fe 



= V M V ( 2fc ) ! ^^(fc-r) r(2r)r(2(fc-r + 2))2 2 - 2 ( fc + 2 ) 

^ o (2fc!) T(fc + 3) ^ (2r)!(2(fc - r))! P r(r)r(fc - r + 2) 

JL V r £ \ 2 * V (2(fc-r) + 3)! 

2 2 ^T(fc + 3) V2J ^ P r!(2(fc-r))!(fe-r + l)! 

^ ~ (-1)^ , £ x 2fe * 2 2(fc _ r) (2(fc-r) + 3)(2(fc-r) + 2)(2(A:-r) + l) 

22Ajr(fc + 3)b; ^ P r!(fc-r + l)! 

tt ~ J-l)*_/ 5 N2 fc * 2 _ r) (2(fc-r) + 3)(2(fc-r) + l) 

2^-r(fc + 3)V27 ^ 1 r\(k-r)\ 



(=o r(fc + 3)V 2 y ^ H(*-r)I 

(_!)* /zx2fc^ r) 4(fc-r) 2 + 8(fc-r)+3 



r!(/c — r)! 



= v (-l) fc P V K V a 2r fl 2 (*- r ) V gj "' 2 
2^r(fc + 3)V2/ ^ ^ Z-j r !(A:-r-j)! 

where a ,2 = 3,ai i2 = 12, 02,2 = 4. Therefore, we obtain that 
exp{z,z(a! cos 9 + (3 sin 9)} sin 4 9d9 

i v ( ~ 1)fc f z -) 2k y ^1 yi( k ~ A cr&k-j-r) 

*y±^( z -) 2k y^^^+^ 
*y a - 2 y ( ' 1)fc ( z -) 2k V— {jtfTp?*-* 
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(-l) fe+J /Z\2(fc+J) 1 



2 oo 



fc! 



'21 



3=0 



_>/,- 



fe=0 



= 2tt J 3 J 2 (zVa 2 +/3 2 ) ^ J 3 {z^a 2 + f3 2 ) + - £ J 4 (z V« 2 + /? 2 ) 



(4.4) 



By means of the same arguments we have that 

e «(a cos +/ 3 S inS) ^6 QdQ 



-Y 

92 



(_1)* (z s2k 



2 2 ^ r(fc + n+l) V2 



a j,3 



r=0 



2tt< 



15 



45/3 2 



r\(k — r — j)\ 



(z^a 2 + /3 2 ) 3 z 2 {^fu 2 + /3 2 ) 4 z(^a 2 +/3 2 ) 5 



P 6 



Mz^/a 2 + P 2 ) 



(4.5) 



where a , 3 = 15, ax,3 = 90, a 2 ,3 = 60, a 3i 3 
For !/ = n, with n > 0, one has that 



z(a cos 9+13 sinS) gm 2n 



^ (2fc!) ^ 

k=Q y ' r=0 



,7T/2 

rg2(*-r) / ( 
V 



cos^^sin 2 ^"^") 



,2fc V- f 2k \ n 2r om-r) ^ + Vg^ ~ r + 71 + 1/2) 



(2fc!) 



T(/c + n + l) 



= *rE 



M (-1)* z 2 ^ 

(2fc!) r(fc + n + l)^(2r)!(2(fc-r)) 



A:=0 



* (2fc)! 2(fc _ r) r(2r)r(2(fc^r + n))2 2 - 2 ( fc +") 

f9,rW9rfc - rV>! P T(r)r(fc-r + n) 



(-l) fe /z 



2 2 "~ 2 ^r(/c + n + l) V2 

* y (-i) fc ^ 

2 2«-2 r(/c + n 

fe=0 v 

7T (~l) fe r .x2A- 

On-l / i 



V ^^(fc-r) (2(fc-r)+2n-l)! 
^ P r!(2(Jfe-r))!(*-r + n-l)! 



r=0 
fc 



fc=0 



fff'V ar fl 2(fc-r) ( 2 ( fc - + 2n - - r) + 2rc - 2) ■ ■ ■ (2(fc - r) + 1) 

V 2 7 ^ P r!(fr-r + n-l)l 

r(fc + n + l) V2/ ^ p 



r=0 



2r R 2(k-r) (gft - *Q + 2n - l)(2(fc - r) + 2n - 3) ■ ■ ■ (2(fc - r) + 1) 

r!(fc — r)! 



Therefore, by taking into account similar manipulations appearing in ( 4.3 1 , (4.4) and (4.51, 
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we can infer that the following equality holds 

(2(fc -r)+2n- l)(2(fc - r) + 2n - 3) • • • (2(fc - r) + 1) _ ^ flj> 

r!(fc — r)! r\(k — r — j)\ 

where a^s are positive constants. We are not able to explicit a^ n for all values of n and j 
because for the above constants seem that a recursive rule does not hold. Nevertheless, it is easy 
to see that 

<*>,„ = (2n - l)(2n - 3) ■ ■ ■ 3 ■ 1, a„,„ = 2" 
Furthermore, the cases tackled above can be sum up in the following scheme 

n = — > ao.o = 1 

n = 1 -> a ,i = 1, oi,i = 2 

n = 2 -> a 0j2 = 3, a lj2 = 12, a 2 , 2 = 4 (4.6) 
n = 3 -> a 0j3 = 15, ai,3 = 90, a 2 , 3 = 60, a 3 , 3 = 8 

Finally, we get that 

I e iz(aco S e+0 S ine) sin 2n edd 

Jo 

oo / -t \k oik n 

2 n-i r(fc + n+ 1) V2/ ^ P l^ r \{k-r-j)\ 

k=0 y ' r=0 j=0 y J ' 

* V (- 1 ) k ( Z \ 2k y Hn^L X^( k -j\^ro2(k-i-r) 

= — V r^ 2 " V aj ' n/32j (a 2 + B*t-i 

- 2 »-^ 8j ^r(H« + i) (fc-j)! lv P) 

j — U k — j 

= — — y(-l) J aj J 23 - PJ J„+,(zVa 2 +/3 2 ) 

= 2?r y (-l poj.n^ ( z ^T^g2) (4.7) 

^2Jz™-J(V a2 + / 32 ) rl+3 

Important results are also the following ones 

Jo V27rr(/x + v + 1) p 2 
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with Re > —\ and Rev > — \ (see Gradshteyn-Ryzhik, 1980, formula 6.581(3)), 

f M^Ma-x)^ = (i + i) -W^) (4.9) 
J Q x a — x \ /i v ) a 

with Re(n) > 0, Re(v) > (see Gradshteyn and Ryzhik, 1980, formula 6.533(2)), and 

(n 2 - h 2 Y^^b^ 

x^- y J u+1 {ax)JJbx)dx = 1 - / (4.10) 

with a > b, Re(v + 1) > Re(/i) > (see Gradshteyn and Ryzhik, 1980, formula 6.575(1)), 

(sinx) 

o 



(sinx)" cos(6cosa;) J^(asinx)(ia; = i / — ? — i , (4.11) 



for Re v > -1 (see Gradshteyn-Ryzhik, 1980, pag. 743, formula 6.688.(2)) 
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